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In this paper, a model for the three-dimensional Spacecraft COntrol 
Laboraory Experiment (SCOLE) is developed. The objective behind this method 
of modeling is to utilize the basic partial dififemtial equations of motion for this 
distributed parameter system and not to use the modal expansion in developing 
the model. The final model obtained is in terms of a transfer function matrix 
which relates the flexible mast parameters like displacement, slope, shear stress 
etc. to external forces and moments. 

1. INTRODUCTION 


It is widely recognized that the future space exploration would require a 
wide array of very large and flexible spacecrafts with very stringent pointing and 
vibration suppression requirements. Some of these spacerafts would also be 
deployed as an assemblage of a number of flexible members. In order to design 
control systems to meet these requirements, accurate dynamical models of the 
flexible spacrafts would have to obtained. Generally, the basic dynamical equa- 
tions are developed in terms of a system of partial differential equations and one 
common approach is to formaulate solutions of these equations in terms of an 
infinite modal expansions and use this approach for developing control systems. 

In this paper, an attempt is made to work with the basic partial differential 
equations and by using Laplace Transforms and incorporating boundary condition 
relationships an alternate modeling scheme is proposed. This methodolgy is based 
on extensive details documented in reference tl] and is applied to NASA Langley 
Research Center s SCOLE problem [2,3]. 6 y 

2. NOMENCLATURE 

u x (t,z) Displacement at point z in roll bending 
Slope of beam at point z in roll bending 
M x (t,z) Bending moment at point z in roll bending 
Shear stress at point z in roll bending 
f x (t,z) External force per unit length in roll bending 


922 



- 2 - 


I x Moment of inertia of area about neutral axis in roll bending 

B x Damping in roll bending 

u x (t,z) Displacement at point z in pitch bending 

Q y (t,z ) Slope of beam at point z in pitch bending 

M y (t,z) Bending at beam at point z in pitch bending 

o y (t, z ) Shear stress at point z in pich bending 

f y (t,z) External force per unit length in pitch bending 

I y Moment of inertia of area about neutral axis in pitch bending 

B y Damping in pitch bending 

xjj f (t,z) The angular displacement at z of an element dz of the beam 

I p The polar moment of inertia of the cross-section 
p The mass per unit volume 
I p G The torsional stiffness of the beam 
G The shear modulus of the material 
I Modulus of elsticity 
p Mass per unit length 
3. METHODOLOGY 

The partial differential equations governing the roll bending motion are 

du X fl (1) 

nr =e * 

This is obtained from the definition of slope. 


36 * 

dz 


(ED. 


1 ~M X 


(El) 


( 2 ) 


This equation is based on beam theory and here (EI) X and (£/), are considered 
equal and represented by (El) 


dM x (3) 

dz 

This is obtained from the definition of bending moment. 
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= PTT- -f x (t t z) 


This equation is based on Newton’s law of motion and the term f x (t,z ) = -B x — 
if there is no external force per unit length on the beam. ™ 

The corresponding equations for pitch bending are 

du y 

dz = e > ( 5 ) 

ae„ , 


= P^r 


The following equations describe torsional bending 

<hL = _J_ T 

dz I p G 

This is based on torsional flexibility. 


l P a,2 


-My(t,z) 


This equation ^is based on Newton’s Law of motion and the term 
M y(t,z) = -By if there is no external torque per unit length of the beam. 

Defining the state variables as 
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the following equations are obtained 

d<7i 


d<72 1 

dz ~ ( EI ) q 3 


dq 3 
dz 


= <74 


=p(!4)«.-a(w) 


dz 


dr 


d<75 

dz 

d<76 

5z 

d<?7 

dz 


= <76 


(£/) 


qi 


= q 8 


( 12 ) 


dq&_ 

dz 


= P(77>75 -/y(^ z > 

dr 


d<?9 _ _1_ 
dz “ I p G 


q to 


= p/ p (-^y)<79 -M v (f,z) 

dz dr 

These equations can be expressed in the form 




03 ) 


Taking Laplace Transforms of the previous matrix-vector equation, the following 
equation is obtained. 
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(fe ~ ( F 0 + F i s + F 2 s 2 )Q - F \q(0,z) - F 2 [ <7(0, z ) + sq(0,z) ] + U(s,z ) 


dQ. 

dz 


= (F 0 + F is + F 2 s 2 ) Q + U(s,z) 


(14) 


Here Fi is a 10 x 10 null matrix. The matrix 
except the following nonzero elements. 


Fq is a 10 x 10 matrix with zeros 


* 0 ( 21 ) = 1 


F o(3,2) 


-1 

(El) 


F o(4,3) = 1 


F 0 (l,6) 


-1 

(El) 


F o(8,7) = 1 


F o(9,9) 


1 

I P G 


The matrix F 2 is also a 10 x 10 matrix with zero entries except for the following 
elements 


^2(4,1) = P 
^2(8,5) = p 
F 2(10,9) = p/ p 

The equation (14) represents a linear system and this can be solved by using 
the state-transition matrix as J 6 

2 

Q(s,z) = H(s,z - z Q ) Q(s,z 0 ) + J H(s,z - %)U {s,\)d% . (15) 

*o 

The SCOLE model is of finite length; i.e. 0 <z<L. Then, at z = L, 

Q(s,L) = H (s,L) Q(s, 0) + J H(s,L - fyU (s&dt, . (16) 
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To determine Q(s,L) and (2(5,0) which in turn would allow determination of 
0(s z) for any t between 0~and L, ten terminal relations must be specihed. these 
terminal relations can be expressed in the form of ten ordinaiy differential equa- 
tions in the following vector-matrix form using linear differential operators. 


M (D) q(t, 0) + N (D) q(t,L) = Fit ) 


(17) 


Here Fit) represents external control forcing functions ( the physical inputs like 
forces'on the flexible beam or shuttle moments etc. or linear combinations thereof 
). Taking Laplace Transforms on both the sides of equation (17), the following 
vector-matrix equation is obtained. 


M is) (2(5,0) + Nis ) QiS,L ) = Fis) 


(18) 


The vector F(s) is the sum of the Laplace Transforms of FQ) and any initial con- 
dition terms. The termination is said to bo homogeneous iFF( 5) - 0. 

It is important to note here that each of the differential equations given in 
(17) in case of a distributed parameter system may involve quantities at both ends 
of the distributed system. This represents termination characteristics or a feed- 
back system where quantities at one end are made to depend on quantities at t e 
other. The boundary conditions for SCOLE model without any external forcing 
functions are given as follows. 

At the shuttle end where z = 0, 


d u x 

iEI ) a 2 (0,0 - o 

Moment 

(19a) 

d It' y 

(EI)—f (0,/) - 0 

Moment 

(19b) 

d u x 

(ED 3 (0.1) -0 

ShearForce 

(19c) 

d « v 

iEI )-r- JL io,t) - 0 
oz 3 

ShearForce 

(19d) 

(7pG)~jp (0,0 = 0 

Torquq 

(19e) 

At the reflector end where z = 

L, the corresponding boundary conditions are 

5 

; cl; 
^ sf 

II 

O 


(20a) 
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3 m v 

(El)— f ( L,t ) = 0 

3 u x 

(ED— f (L,t) = 0 

3 My 

(£7)-^ (L,0 = 0 
( /pG )^ = 0 


(20b) 

(20c) 

(20d) 

(20e) 


Hnn r n^ der i?i° 0 ? ain a complete representation of the system governed by equa- 
tions (16) and (18), we substitute (16) into (18) with z 0 = 0 and z = L\ * ^ 


[ M{s) +N(s ) H(s,L) ] Q(s,0) = V^y) 


where, 


L 


V(s) = F(s) - N(s)j H(s,L - fyUis&dt, 
0 


( 21 ) 


( 22 ) 


ST I ? / d ? a /. the relationship given in (17) is independent and as a result 
[ M(s) + N(s)H(s) ] is of rank 10 for SCOLE and so has an inverse. Hence we 
can write 9 


Q(-y.O) = [ M(s) + N(s)H(s,L ) ] 1 V%y) (23) 


Thus, 


Q(j,z) = H(s,z) [ M{s) + N(s)H(s,L) ] K(j) + J H(s,z-^)U(s^)d^ (24) 

0 

If there are no distributed forcing terms and initial conditions, equation (24) can 
be written as n v ' 


Q(s,z) = H(s,z ) [ M(s ) + N(s)H(s,L ) ] _1 F(s) 


(25) 


where F(s) is the Laplace transform of the external forcing terms in (16). Hence 
the mamx of transfer functions from F(s) to Q(s,z) is given by 


G(s,z) = H(s,z ) [ M(s) + N(s)H(s,L ] 


-l 


(26) 
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